Critical behavior and driven Monte Carlo dynamics of the XY spin glass in the phase 

representation 



o 
o 

(N 
ON 



o 
o 



> 

NO 
(N 
(N 
O 

o 

-I— > 
c3 



o 
o 



X 



Enzo Granato 
Laboratorio Associado de Sensores e Materials, 
Instituto Nacional de Pesquisas Espaciais, 
12245-970 Sao Jose dos Campos, SP Brazil 

A driven Monte Carlo dynamics is introduced to study resistivity scaling in XY-type models in 
the phase representation. The method is used to study the phase transition of the three-dimensional 
XY spin glass with a Gaussian coupling distribution. We find a phase-coherence transition at finite 
temperature in good agreement with recent equilibrium Monte Carlo simulations which shows a 
single (spin and chiral) glass transition. Estimates of the static and dynamic critical exponents 
indicate that the critical behavior is in the same universality class as the the model with a bimodal 
coupling distribution. Relevance of these results for 7r-junction superconductors is also discussed. 
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Phase transition in the three-dimensional (3d) XY spin 
glass with short-range interactions has been a challeng- 
ing problem. The additional chiral order parameter 1 , 
with Ising-like symmetry, arising from frustration effects 
competes with spin glass ordering due to freezing of the 
two component spins. Earlier work showed that while 
the spin glass transition temperature vanishes, the chi- 
rality orders at a finite temperature^. Evidence of such 
behavior has been provided both for models where the 
coupling between XY spins has a bimodal distribution ( 
+1 and —1 with equal probability) and for models with a 
Gaussian coupling distribution. However, recent numer- 
ical work based on equilibrium and dynamic simulations 
has questioned this decoupled scenario with some con- 
flicting results. 

Like other XY type models^, the XY spin glass can 
also be regarded as a model for phase-coherence in su- 
perconductors. In particular, it is currently being used 
as a model for granular superconductors containing ir 
junctions^SiSiS, as in high-T c superconductor materi- 
als with d-wave symmetry^. In this case, a tt junction 
corresponds to an antiferromagnetic coupling between 
the XY spins while the orientational angle of the spins 
represents the phase of the local superconductor order 
parameter. This superconductor analog allows the use of 
the electrical resistivity as a very useful dynamic quantity 
that can characterize spin-glass ordering in the XY spin 
glass. In fact, the linear resistivity is a direct measure 
of phase stiffness and therefore phase coherence in the 
superconductor, which is equivalent to long-range order 
in the spin variables. Thus, the resistivity behavior can 
be used to study numerically the phase transition in the 
XY spin glass. 

Earlier Monte Carlo (MC) simulations of resistivity be- 
havior of the XY spin glass with a bimodal coupling dis- 
tribution in the vortex representation's, showed evidence 
of a resistive transition at finite temperature. However, 
this result was interpreted as an indication of the chiral- 
glass transition. With a different interpretation^, it was 
argued that this resistive transition should instead be at- 



tributed to the spin-glass ordering. The possibility of 
spin-glass ordering at finite temperatures was also sup- 
ported by calculations of the spin stiffness exponent in 
the ground stat o 11 ! 12 , showing that the lower-critical di- 
mension is below 3, which implies that a phase-coherence 
transition at finite-temperature is possible in three di- 
mensions. Later on, resistivity calculations in the phase 
representation using Langevin dynamica&S confirmed the 
occurrence of a phase-coherence transition at finite tem- 
perature and suggested a single transition scenario where 
chiral and phase variables order simultaneously. 

Recent equilibrium MC simulations by Kawamura 
and Li for the XY spin glass with a bimodal coupling 
distribution^, provided an improved estimate of the 
chiral-glass transition temperature, which turned out to 
be consistent with the transition temperature as obtained 
from the resistivity scaling-, but still conclude for the 
absence of phase coherence within a spin-chirality de- 
coupling scenario. However, very recently, MC calcu- 
lations of the chiral and spin correlation lengths by Lee 
and Young for the model with Gaussian couplings^ show 
that there is indeed a transition at finite temperature at 
which both spin and chirality order, supporting the single 
transition scenario. These conflicting results raise the im- 
mediate question of universality class of the phase transi- 
tion in the two models. In addition, since the resistivity 
scaling approach indicates a phase-coherence transition 
for the bimodal coupling model at the chiral-glass tran- 
sition, it should also be of particular interest to find out 
if the gives consistent results for the Gaussian coupling 
model. 

In this work, we study the critical behavior of the 3d 
XY spin glass with Gaussian couplings by a scaling anal- 
ysis of the nonlinear resistivity obtained by MC methods. 
First, we introduce a driven MC dynamics for XY-type 
models in the phase representation. The alternative MC 
dynamics in the vortex representation**, which has been 
used for the model with bimodal couplings, is not use- 
ful for the Gaussian model because in this case the cou- 
pling magnitude is not uniform. The main advantage of 
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the driven MC method in the phase representation com- 
pared with standard Langevin dynamics simulations 8 9 
is that much longer time scales can be accessed, insuring 
that the long-time behavior is probed at the lowest tem- 
peratures and current densities. In addition, the driven 
MC dynamics is particularly useful in models for which 
Langevin phase dynamics or MC vortex dynamics are not 
available^. 

From a scaling analysis of the resistivity data for the 3d 
XY spin glass, we find a phase-coherence transition at fi- 
nite temperature and the corresponding thermal and dy- 
namic critical exponents are determined. The results are 
in good agreement with recent equilibrium MC simula- 
tions for the model with Gaussian couplings 14 . Moreover, 
comparing the static and dynamic critical exponents ob- 
tained for the Gaussian model with recent results for 
the bimodal coupling model using the same analysis^ 
suggests a common universality class. Altogether, these 
results strongly support the single transition scenario 8 
where chiral and phase variables order simultaneously in 
both models. 

We consider the XY-spin glass driven by an external 
perturbation, described by the Hamiltonian 

H = -J2 -h cos(0< - fly) - jYJfii ~ 6 i+x) (!) 

<ij> i 

The first term gives the nearest-neighbor coupling energy, 
where Jij are quenched random variables taken from the 
Gaussian distribution with zero mean and unit variance 

P(J ij ) = -^=exp(-J? j /2). (2) 

The second term in Eq. represents the effects of an 
external perturbation, applied in the x-direction for con- 
venience, coupled to the bond phase difference &i — 0i +x . 
When regarded as a model for granular superconductors 
with 7r junctions*!* 1 ^, the random distributed negative 
signs of the bond variable Jij correspond to Josephson 
junctions with a phase shift of n and the perturbation 
is equivalent to a driving current density J applied to 
the superconductor. When J ^ 0, the total energy is 
unbounded and the system is out of equilibrium. The 
lower energy minima occur at increasing phase differ- 
ences 0i — 9i+ x as a function of time, leading to a net 
phase slippage rate proportional to < d(9i — 9i +x )/dt >, 
which can be taken as a measure of the voltage V (in 
arbitrary units) in a model of superconductors. 

To study the nonequilibrium behavior generated by the 
driving current density J in Eq. ^ we use a driven MC 
dynamics method. The time dependence is obtained by 
identifying the MC time as the real time t and we set the 
unit of time dt = 1, corresponding to a complete MC pass 
through the lattice. Periodic (fluctuating twist) bound- 
ary conditions are usecU^ in cubic systems of linear size L. 
This boundary condition adds new dynamical variables, 
u a (a = x,y and z), corresponding to a uniform phase 
twist between nearest-neighbor sites along the principal 



axis directions x,y and z. A MC step consists of an at- 
tempt to change the local phase 9i and the phase twists 
u a by fixed amounts, using the Metropolis algorithm. If 
the change in energy is AH, the trial move is accepted 
with probability min{l, exp(— AH/kT)}. The external 
current density J in Eq. ^ biases these changes, leading 
to a net voltage (phase slippage rate) across the system, 
given by 

1 d L 

V = J2 J t E " hd* ~ U * L ^ ( 3 ) 

Using this procedure, the electric field E — V/L and 
nonlinear resistivity p = E / J can be obtained as a func- 
tion of the driving current density J for different tem- 
peratures T. The main advantage of this MC method 
compared with standard Langevin dynamics 7 ^ for mod- 
els in the phase representation is that in principle much 
longer time scales can be accessed. In the latter method, 
the maximum time step is limited by numerical instabili- 
ties when integrating the Langevin differential equations 
whereas in the former, trial moves are possible which 
would correspond to very large time steps in the numer- 
ical integration. 

For the simulations of resistivity behavior described 
below, first MC calculations are performed with J = 
(zero current bias) to obtain the equilibrium state which 
is then used as initial state for the driven MC dynamics 
with J 0. Extensive calculations were performed, using 
typically 10 7 MC steps for the equilibration and driven 
MC dynamics with 10 to 20 different realizations of the 
Jij disorder distribution for low values of J in system 
sizes ranging from L = 4 to L = 12. 

To extract the critical behavior from the numerical re- 
sults of the nonlinear resistivity p we need a scaling the- 
ory for the resistive behavior near a second-order phase 
transition. A detailed scaling theory has been described 
in the context of the current-voltage characteristics of 
vortex-glass models^ but it can be directly applied here. 
If a phase-coherence transition occurs at nonzero tem- 
perature T c , then measurable quantities scale with the 
diverging correlation length £ oc \T/T C — and relax- 
ation time t a f , where v and z are the thermal and 
dynamical critical exponents, respectively. The nonlinear 
resistivity p should then satisfy the scaling form 3 

TK^^fty) (4) 

where g(x) is a scaling function. The + and — signs 
correspond to T > T c and T < T c , respectively. If the 
numerical data satisfy such scaling form for different tem- 
peratures and driving currents, then the critical temper- 
ature and critical exponents of the underlying equilib- 
rium transition at J — can be estimated from the best 
data collapse. However, for a reliable estimate, the data 
should also satisfy the expected finite-size behavior in 
smaller system sizes. Finite-size effects are particularly 
important sufficiently close to T c when the correlation 
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FIG. 1: Nonlinear resistivity p for different temperatures T, 
for a system size L = 12. 



FIG. 2: Scaling plot of the data from Fig. 1 near the transi- 
tion, with £ oc \T/T C - l\~ v . 



length £ approaches the system size L. In particular, 
at T c , the correlation length will be cut off by the sys- 
tem size in any finite system and the nonlinear resistivity 
should then satisfy a scaling form as in Eq. 0]with £ = L, 

JT2 

T cP L^=g{—) (5) 

Away from T c , the scaling function in Eq. [3|will also de- 
pend on the dimensionless ratio^iS L/£ as g( J£ 2 /T,L/£). 
To simplify the analysis, we consider resistivity data at 
current densities such that J£ 2 /T = is constant. Then, 
the scaling form depends only on a single variable and 
the resistivity should satisfy the finite-size scaling form 

TpL*^ = g{L 1 ' v {T/T c -l)) (6) 

We now describe the results of the scaling analysis of 
the resistivity behavior obtained from the driven MC dy- 
namics. The nonlinear resistivity p as a function of cur- 
rent density and temperature is shown in Figs. 1 for a 
large system size L = 12. The behavior is consistent 
with a phase-coherence transition at an apparent critical 
temperature in the range T c ~ 0.25 — 0.35. At higher 
temperatures, the linear resistivity p^ = Voiij^qE/J is 
finite while at lower temperatures, it extrapolates to very 
low values. The actual limiting values of the resistivity 
at low currents can not be determined accurately with 
the available computer time. However, if one assumes 
a continuous equilibrium transition at a finite tempera- 
ture with J = 0, then the nonlinear resistivity should 
satisfy the scaling form of Eqs. 0]and the critical tem- 
perature and exponents can then be obtained from the 
data collapse. Such scaling plot is shown in Fig. 2, ob- 
tained by adjusting the unknown parameters. We now 
show that theses estimates, using the largest system size, 
are reliable by verifying that they satisfy the expected 
finite-size behavior using smaller system sizes. In fact, 
as shown in Fig. 3, the nonlinear resistivity satisfy the 
finite-size scaling form of Eq. [SJfor different system sizes 



at the estimated T c — 0.335. Also, away from T c , the 
nonlinear resistivity calculated at different temperatures 
and system sizes satisfy the scaling form of Eq. as 
shown in Fig. 4. From this scaling analysis we obtain 
the estimate of critical temperature and critical expo- 
nents T c = 0.335(15), z = 4.5(3) and v = 1.2(2). It 
should be emphasized here that this result is obtained 
by requiring that T Cl z and v satisfy, not only the scal- 
ing form of Eq. 0] for the large systemii, but at same 
time also the finite-size scaling forms of Eqs. and for 
smaller systems. 

Our estimate of the phase-coherence critical tempera- 
ture and exponent for the model with Gaussian couplings 
is particularly interesting in view of the recent equilib- 
rium MC simulations 14 which show evidence of spin glass 
transition at finite temperature and that chirality and 
spin variables order simultaneously. Our results provide 
further support for the estimate of the spin-glass tran- 
sition at T c — 0.34(2), as first reported in that work. 
Moreover, our estimate of the static and dynamic critical 
exponents agree with those obtained for the model with 
bimodal coupling distribution using similar analysis^ ( 
v = 1.2(2), z = 4.4(3) ), suggesting a common static and 
dynamic universality class. If indeed the critical behav- 
ior is the same for both models then a single transition 
should also be observed in the latter model, despite the 
conclusions from other MC simulations 1 ^ that find a spin- 
chirality decoupling scenario. 

Our results for the resistive transition in the Gaus- 
sian coupling model are also relevant for understand- 
ing the behavior of granular superconductors with 
7r junctions^SiSiS. In particular, some measurements 
in high-T c superconductor materials 1 ™, and numerical 
simulations 1 ^, showing a power-law behavior for the non- 
linear contribution to the resistivity near the onset of the 
paramagnetic Meissner effect, have been interpreted as 
resulting from a chiral-glass transition with no phase co- 
herence, based on earlier results for the XY-spin glass 
model with a bimodal coupling distribution. In view of 
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FIG. 3: Finite-size scaling plot of the nonlinear resistivity at 
T c = 0.335. 
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FIG. 4: Finite-size scaling plot near T c using current densities 
such that J£ 2 /T = 1, a constant value. 



the results for resistivity scaling, an alternative interpre- 
tation is possible^ where the observed behavior is a conse- 
quence of the underlying phase-coherence transition and 
the power-law exponent a is determined by the dynamic 
critical exponent z as a = (5 — z)/2. In this regard, the 
results for the Gaussian model considered here have two 
important implications. First, this model is a more re- 
alistic description of the superconductor material since 
in a granular sample the Josephson couplings between 
grains will have random magnitudes, as well as random 
signs due to the it phase shifts. Secondly, our results 
show that a phase-coherence transition takes place at fi- 
nite temperature and the common universality class we 
find for the resistive transition in these models shows that 
the proposed^ superconducting chiral transition and the 
power-law exponent a, do not depend on details of the 
coupling distribution. 

In summary, we have introduced a driven MC dynam- 
ics method to determine the resistivity behavior of XY 
type models in the phase representation. The method 
is used to study resistivity scaling and the phase transi- 
tion in the 3d XY spin glass. From the scaling analysis 
we find clear evidence of a phase-coherence transition at 
finite temperature. The critical temperature and expo- 
nents are in good agreement with recent equilibrium MC 
simulations with a Gaussian coupling distributions^ and 
suggest that the critical behavior is in the same univer- 
sality class as the the model with a bimodal coupling 
distribution*. For 7r-junction superconductors, the re- 
sults demonstrate that the superconducting chiral-glass 
transition and the numerical value of power-law expo- 
nent a, are not sensitive to the details of the coupling 
distribution. 

This work was supported by FAPESP (grant no. 
03/00541-0). 
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